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Abstract 

Three magnon bound state problem is studied within Bethe Ansatz for an exactly 
rung-dimerized spin ladder. It is shown that contrary to the general three-magnon 
problem the solvability is less sensitive to non integrability. All obtained wave 
functions are presented in explicit forms. 

1 Introduction 

During the last century a wide number of methods were developed for analysis of quantum 
integrable systems [H |2l |3]. However almost all physically interesting models are non 
integrable. It is always assumed that in this case all the analytical methods (which may 
be successfully applied for integrable models) should fail. In the present paper we show 
that the situation is not so extremely hopeless. Treating the general model of exactly 
rung-dimerized spin ladder [H |3] we obtain general conditions for solvability of three- 
magnon bound state problem within Bethe Ansatz. Surprisingly it turned out that these 
conditions are essentially weaker than the corresponding integrability ones. 

Bethe Ansatz approach to three-magnon states in an exactly rung-dimerized spin lad- 
der |1] was previously studied by the author [5]. It was shown that the general problem 
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is completely solvable only in five special integrable cases. Solvability of a more special 
bound state problem at the first glance also implies integrability. Indeed a bound state 
may be obtained [H El 12] from scattering ones by analytic continuation of the corre- 
sponding wave numbers. However one may choose another way suggesting Bethe Ansatz 
for bound states as a separate problem when only bounded exponents should be leaved. 
Surprisingly this approach is less sensitive to non integrability. 



2 Spin ladder Hamiltonian 

Spin ladder Hamiltonian acts on an infinite tensor product of spaces related to the ladder 
rungs (numerated here by an index n) 

n = l[mn. (1) 

n 

For each n one has hn = C^ = C^(S) where the two C^-factors are representation spaces 
for two triples of 5 = 1/2 spin operators Sj^n (j = 1, 2) associated with the n-th rung 



Si,nSi,n — ^^afecSj „, [S" „, 83 „] — 0, a,b,c — 1, 2, 3. (2) 



With respect to action of the total rung-spin operator 

Sn = Si_„ + S2,n- (3) 

each space hn may be decomposed into three-dimensional rung-triplet and one- dimensional 
rung-singlet subspaces 

K = hi® (4) 

In the present paper we use the following representation for Hamiltonian of an exactly 
rung-dimerized spin ladder [5] 

n 

+ J4S ) • (5) 

Here 

o„ = 

2 



Qn = (6) 



is rung-triplet projector 

Qn\hf^ = 0, Qn\hi = 1- (7) 
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The two Hermitian conjugated su(2)-covariant operator triples [5] 

*n = ^(Sl,„ - S2,„) - i[Si,„ X S2,„,], 

*n = ^(Sl,.-S2,„)+«[Si,, X S2,„], (8) 

additionally to 

[S^, = tsabc^l [SI = tSabA a,b,c= 1, 2, 3, (9) 
satisfy relations 

S:*^ = 0, nS'n = 0. (10) 

Both (9) and (10) may be proved from (2). According to (6), (9) and (10) 

So and may be treated as (neither Bose no Fermi) creation-annihilation operators 
for rung-triplets. The subspace is generated by a single vector |0)„ for which 

*„|0)„ = 0, S„|0), = 0, Qn|0)„ = 0. (12) 

We shall utilize the following two basises of the space /i^. The covariant basis 

|1): = *^|0)„, a = 1,2, 3, (13) 

for which according to (9) 

S^|l)^ = ^e,,e|l);;, (14) 

will be used for representation of three- magnon states with total spin and 1. The 
standard S^-eigenbasis 



s^|i)^ = 0, s^m = V2\iZ s^\iZ = V2\i)t 

sl\i)^ = ±|i)±, s^|i)° = o, (15) 



will be used for representation of three-magnon states with total spin 2 and 3. According 
to (9) one can suggest the following representation 

|i)^ = T*^|o)n, |i)° = |i)^, (16) 

where 

*±^i=(*i±z*^j^ ^^^^[^i±t^l). (17) 



The formulas 

*^|l)^ = 5,b|0)„, a, 6, = 1,2, 3, (18) 

may be obtained from (13) by Hermitian conjugation. Analogously for the operators 
one have 

*±ri)^ = ±io)„, *^ri)^ = *^ri)° = 0. (i9) 

Eqs. (7) and (11) result in the commutation relation 

[H,Q]=0, (20) 

where 

q = J2Q^^ (21) 

n 

is (according to Eq. (7)) the number operator for rung-triplets. As it follows from (20) 
an infinite tensor product of rung-singlets 

\0)r-d = ll^\0)n, (22) 

n 

is an eigenvector for H. It was already mentioned [H [5] that for rather big Ji > it is 
the ground state of the system. In this (exactly rung-dimerized) case the physical Hilbert 
space is an infinite direct sum of multi-magnon sectors 

CXI 

-^Phys ^ ym^ g|^^ ^ ^_ ^23) 

m=0 

Interpretation of the Hamiltonian (5) in this case is clear. The first term describes 
the chemical potential of an excited rung-triplet. The second one corresponds to rung- 
triplets kinetic energy. We shall imply that J2 7^ 0. The last three terms describe a 
spin-dependent interaction between two neighboring rung-triplets. 



3 General properties of three-magnon wave functions 

We shall use the notation 15, k) for a three-magnon state with total spin S. At 5 > it 
will be supplied by an upper index. For representation of a S* = 1 state |1,A;)" will be 
used the basis (13) (so in this case a = 1, 2, 3). At S* = 2, 3 using the basis (16) we shall 



4 



represent only the vectors \S, k)^ with = S. Namely we shall treat the following states 




\0,k) = eatc Yl e^'^™+"+^')/36o(fc,T^-m,p-n)...|l);:,...|l)^..|l)^ 

m<n<p 

\l,kr = Yl e^'^™+"+^)/3[6«(A;,n-m,p-r^)...|l);:,...|l)^..|l)J. 



m<n<p 

'P) (u ^ ^ ^ ^ rr,\ h\a 

p ' 



+ bY>ik,n-m,p-n)...\l)l...\ir„...\l) 
+ b?{k,n-m,p-n)...\l)l...\l)l..\l);... 
\2,k)' = Yl e^'^-^-^'^/'[bi'\k,n-m,p-n)...\l)i...{\l^ 

m<n<p 
d2) 



+ bf{k, n - m,p - n) . . . . . . - . . . . . . |i) 



V 



|3,fc)3 = Y e^'^"^+"+^)/%(fc,n-m,p-n)...|i)+...|i)+...|i)+... (24) 



m<n<p 



All the (reduced to center mass) wave functions bj{k,m,n) (j = 0,1,2,3) have a 
physical sense only at m, n > 0. However within the Bethe Ansatz approach [3] they 
should be continued into the two unphysical boundary regions (m = 0, n > 0) and (m > 
0,n = 0). At m, 77, > 1 the corresponding Schrodinger equations have an identical form 


6Jibj{k, m, n) + J2[e~^''^%j{k, m + 1, n) + e^^/%j{k, m — l,n) + e~'^^/^bj{k, m — 1, n + 1) 
+e*^/^6j(/c, m + 1, n — 1) + e~''^/%j{k, m, n — 1) + e'-^/%,j{k, m, n + 1)] = Ebj{k, m, n)(25) 

for all j = 0, 1, 2, 3. Requiring correctness of (25) in the boundary regions (m = 1, n > 0) 
and (m > 0, n = 1) one immediately obtains the following systems of Bethe conditions [5] 

e'''/%{k,0,n)+e-'''/^bo{k,0,n+l) = 2Ai6o(fc, 1, n), 

e-'''/\{k,m,0) +e'''^\{k,m + l,0) = 2Aibo{k,m,l), (26) 

''^/^foW (A;, 0, n) + e-^'^/^^S') q, n + 1) = ( A2 + A,)b^^^ {k, 1, n) + (A2 - A,)bf^ {k, 1, n), 

"'/^f^ (k, 0, n) + e-^^^^f^ (k, 0, n + 1) = (A2 + A,)bf^ {k, 1, n) + (A2 - A,)b^^^ {k, 1, n), 

'^/%^^\k,m,0) + e'^/%f\k,m + 1,0) = {A2 + A^)bf\k,m,l) + {A2 - Ai)bf\k,m,l), 

"'/¥^\k, m, 0) + e"'/¥^\k, m + 1, 0) = (A2 + Ai)6f ^ {k, m, 1) + (A2 - Ai)6f ^ (fc, m, 1), 



e 
e 
e~ 
e~ 
e 



'"/^^^ {k, 0, n) + e-'^/'b? (A;, 0, n + 1) = ^ ( Aq - A2) [^i'^ (A;, 1, n) + bf^ {k, 1, n)] 

+ 2Aobf\k,l,n), 



5 



e-'^'^h^^^ {k, m, 0) + e'^/^^^^ {k,m + 1,0) = ^ ( Aq - A2) [bf^ {k, m, 1) + bf^ {k, m, 1)] 

+ 2Aob^l\k,mA), (27) 

e'''/¥^^ {k, 0, n) + e-'^Hp {k, 0,n+l) = 2Ai4^^ {k, 1, n) + (A2 - Ai)6^^^ (y^, 1, n), 

e'''/¥^\k,0,n) + e-'''/¥2\k,0,n+l) = 2A26^'^(A;, 1, n), 

e-'^H^^\k, m, 0) + ^^H^^\k, m + 1, 0) = 2Aib^^\k, m, 1) + (A2 - Ai)bf\k, m, 1), 

e-^'^/^fcf (A;,m,0)+e**^/=^6f (A;,m + 1,0) = 2A26^'^(A;, m, 1), (28) 

where 

A„.A^1^, A.^^i^. A.^M_^. p9) 

Z(72 

The system related to b'i{k,m,n) has the form identical identical to Eq. (26) however 
with Ai replaced on A2. 

The systems (26)- (28) are invariant under the set of corresponding duality transfor- 
mations 

'D{bo^3{k,m,n)) = bo^3{k,n,m), V{b^\k,m,n)) ^b^^~^\k,n,m), 

V{b^^\k,m,n)) = b^^'^\k,n,m). (30) 

Since in the present paper we are interesting only in bound states we imply a normal- 
ization condition 

J] |6,(A;,m,n)|2<oo, j = 0,l,2,3. (31) 

m,n>0 

One may ask a question in what interval lies k. On the one hand there should be 

< A; < 27r, (32) 

on the other a substitution k + 27r evidently changes the exponential factors in (24). In 
order to clarify the situation wc notice that this substitution together with simultaneous 
multiplication of the wave function on the factor e^'^*^"*"'*)/^ retains invariant all the states 
in (24) 



4 Bound states at S = (S = 3) 

The following substitution 

bo(k,m,n) ^ z'^z'l, (33) 



solves Eq. (25) resulting in a dispersion 

E{k) = 6 Ji + J2 y^/' (^ + - + z^+ e-^'/' f^ + ^ + zA]. (34) 

\Z- / \Z+ Z- /J 

A substitution of (33) into Eq. (26) gives 

2Aiz± = e±^*^/3 + e'^'''/h^. (35) 

Prom Eq. (35) readily follows 

2Ai(z_ - z+) = e"'/%z+ - Z-). (36) 

In the general case 

2|Ai| ^ 1, (37) 
Eq. (36) results in an autoduality (with respect to (30)) condition 

z- — z+ = z, (38) 

and the system (35) reduces to an equation 

Gi = 0, (39) 

where 

Gj = 2AjZ - ^^'^ - e-'^'^z. (40) 

Eqs. (39), (40) result in 



A substitution of (38) and (41) into (34) gives 

2 7 

E{k) - 6-^1 + 4^/ ^ (8^1 + cos k) . (42) 
Condition (31) is satisfied only if \z\ < 1 or according to (41) 

Aicos/c < (4A^ -3)A2. (43) 

The non-autodual states at 2|Ai| = 1 will be studied in a separate paper. Here we 
shall treat only autodual wave functions. 

As it was mentioned above all calculations for 5" = 3 may be performed in the similar 
manner after the replacement Ai on A2. 
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5 Bound states at S = 1 



Treating an autodual substitution 

bi\k, m, n) 



one gets from Eq. (27) 



V 



Gi + Gs Gs-Gi 
-Gi Gi + Gs 

Gn — Go Gn — Go 3Gn 



5(2) 
5(3) 



(44) 
(45) 



(46) 



Since determinant of the matrix in Eq. (46) is proportional to G0G1G2 the system 
(46) is solvable in three cases 



G,=0, J = 0,1,2. 



(47) 



The corresponding solutions in general do not satisfy Eq. (45). However this problem 
may be solved if the system (46) has two different solutions. Taking into account that 
the latter is possible only if Aj — Ai for j ^ I we conclude that a three-magnon bound 
state related to the wave function exists at three cases 



(48) 
(49) 
(50) 



Ao 


= Ai, 


B = 


[1,-1,1], 


Ao 


= A2, 


B = 


[1,1,1], 


Ai 


= A2, 


B = 


[1,-4, 1], 



where each B 
(46). 



[5(^\ 5(^\ 5(^)] represents the corresponding autodual solution of Eq. 



6 Bound states at S = 2 

An autodual substitution 

b?{k,m,n) = C^^^z"^^", (51) 

C(2) = C^^\ (52) 
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results in 

= 0. (53) 




The system (53) is solvable under one of the conditions (47) taken for j = 1, 2. However 
autodual solutions [1, 1] and [i, —i] exist only at A2 = Ai. 



Ao 


= Ai 


= A2, 






Ao 


= A2 


= ±1, 


Ai 


= 0, 


Ao 


= A2 


= 0, 


Ai = 


±1, 


Ai 


= A2 


^±'-. 
2 


• Ao 


= 0, 


Ai 


= A2 


= 0, 


Ao = 


±1 

2 



7 Summary and discussion 

As it was previously shown in Ref. 5 the system related to Hamiltonian (5) is integrable 
only in five special cases 

(54) 
(55) 
(56) 
(57) 

(58) 

We see that restrictions (48)- (50) are essentially weaker than (54)- (58). So the bound 
states of Bethe form exist even for a rather big number of non integrable models. 

In order to clarify a relevance of conditions (48)-(50) to physical applications we use 
an equivalent representation of the Hamiltonian (5) [5] 

where 

Hn,n+1 — 2^^^'" ' ^^'^ ■ ^2,n+l), -f^n,n+l = ^1," ' Sl,n+1 + S2,n • S2,n+1, 

= Si^n ■ S2,n+1 + S2,n ' Si^„+i, H'!^^^j^'^ = (Si^„ ■ S2,n)(Sl,n+l ■ S2,n+l), 
H^n,n+1 — ' Sl,n+l)(S2,n " S2,n+l)! -^^n/rt+l — (^l.n ' S2,n+l) (82,71 ' Si^„+i). (60) 

A correspondence between the coupling constants of (5) and (59) is given by Eq. (20) of 
Ref. 5. We imply also the condition [U |5] 

Jii-Jdd = ^{Ji-Jd), (61) 
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which guarantees the commutation relation (20). According to Eq. (32) of Ref. 5 

Ao = Ai ^ 12Jz - 8Jd - 5Jh = 0, (62) 
Ao = A2 ^ Ui-Ja = Q, (63) 
Ai=A2 ^ AJd + Jii = Q. (64) 

It is assumed [6] that for physically relevant models there should be 

\Jrr\-, \ Jll\-, \ Jdd\ < \ Jl\-, (65) 

SO Eq. (61) may be satisfied only for highly frustrated ladders with 

Jd ~ Jl. (66) 

In this context the case (62) seems to be physically more reliable. 
The author is very grateful to P. P. Kulish for helpful discussions. 
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